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ABSTRACT 



Context. In the context of the space-based mission CoRoT, devoted to asteroseismology and search for planet transits, 
we analyse the accuracy of fundamental stellar parameters (mass, radius, luminosity) that can be obtained from aster- 
oseismological data. 

Aims. Our work is motivated by the large uncertainties on planetary radius determination of transiting planets which 
are mainly due to uncertainties on the stellar parameters. Our goal is to analyse uncertainties of fundamental stellar 
parameters for a given accuracy of oscillation frequency determination. 

Methods. We generate grids of equilibrium models of stars and compute their pulsation spectra based on a linear nona- 
diabatic stability analysis. Using differents methods of comparison of oscillation mode spectra, we derive uncertainties 
on fundamental stellar parameters and analyse the effect of varying the number of considered modes. 
Results. The limits obtained depend strongly on the adapted method to compare spectra. We find a degeneracy in the 
stellar parameter solutions, up to a few % in mass (from less than 1% to more than 7% depending on the method used 
and the number of considered modes), luminosity (from 2% to more than 10%) or radius (from less than 1% to 3%), 
for a given pulsation spectrum. 
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1. Introduction 

The successful launching of CoRoT produces data for tran- 
siting exoplanets with unprecedented accuracy compared 
to ground-based observations, providing a measure of plan- 
etary mass and radius and thus information on the mean 
density and bulk composition of exoplanets. Unfortunately, 
a remaining source of uncertainty on the planetary parame- 
ters (mass and radius) is due to uncertainties on the stellar 
parameters. 

The precise determination of gravity by spectroscopy 
is difficult, and particularly for F stars. The determina- 
tion of the mass of the star is the main uncertainty on 
the planet's size. This uncertainty on the mass is generally 
about 10% (i.e 3% on the radius and 10% on the density). 
This may be of importance to disantangle a massive planet 
from a brown dwarf, as in the case of the recently discov- 
ered "super- Jupiter" CoRoT-Exo-3b (Leconte et al. 2009). 
It may also prevent distinguishing an ocean planet from an 
earth- like planet (Grasset et al. 2009). 

Our goal is to examine the accuracy on the latter pa- 
rameters , namely mass, luminosity and radius, that could 
be obtained using asteroseismology with the expected ac- 
curacy on oscillation frequencies of CoRoT. We explore the 
space of stellar parameters, in terms of mass, effective tem- 
perature, luminosity, metallicity and mixing length param- 
eter and we analyse the sensitivity of predicted spectrum of 
oscillation frequencies to these parameters. Adopting vari- 
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ous levels of stellar noise and white noise, we analyse the 
frequency uncertainty for oscillation modes of a given am- 
plitude and lifetime. This analysis is based on the per- 
formances of CoRoT in the Asteroseismology and Planet 
Finder channels described in Auvergne et al. (2009). 



2. Method 

2.1. Stars studied 

Among all the transiting systems discovered so far, the large 
difference between the fundamental parameters of XO-3 
given by Johns-Krull et al. (2008) (M=1.41 ± 0.08 M 
; R=2.13 ± 0.21 i? Q ) and by Winn et al. (2008) (M=1.213 
± 0.066 Mq ; R= 1.377 ± 0.083 R Q ) was a strong moti- 
vation to examine how accurately fundamental parameters 
of a star (mass, radius, luminosity) can be constrained by 
asteroseismology. In this preliminary study, we first analyse 
stars with characteristics close to the ones of XO-3 (M ~ 
1.4 Mq , R ~ 1.5 R Q , Teff ~ 6400 K). We focus on modes 
with low order I (0 < I < 3) and radial order n between 5 
and 20. 



2.2. Frequency computation 

The pulsation calculations are performed with a nonradial 
code originally developed by Lee (1985) and based on a lin- 
ear non-adiabatic stability analysis. The equations are lin- 
earised around hydrostatic equilibrium, and eigenfunctions 
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are expressed with spherical harmonics Yi m . The eigenfre- 
quencies are defined by 

a = cr r + i(Ji, (1) 

with ay the oscillation frequency and <7i the damping 
rate (if positive) or growth rate (if negative) (see details 
in Mulet-Marquis et al. 2007, and references therein). The 
system of equations is solved with a Henyey-type relaxation 
method. Stellar structure models are calculated using the 
Livermore opacities (Iglesias, Rogers, 1996). The number 
of gridpoints of the equilibrium models (4000 gridpoints) 
is high enough for the uncertainty on the frequency de- 
termination by the linear code to be less than the accuracy 
expected from CoRoT for a pulsation mode with an infinite 
lifetime and observed for 150 days (i.e 0.1 fiHz). 

2.3. Accuracy of frequency determination 

The photometric performance of CoRoT, described in 
Auvergne et al. (2009), allows for the detection and char- 
acterisation of Sun-like oscillations of bright (6 < V < 9) 
stars in the asteroseismology (AS) channel, and larger am- 
plitude red giant pulsations on fainter (V > 11.5) stars in 
the planet finding (PF) channel. 

The precision to which mode frequencies can be deter- 
mined depends on the mode lifetime T, the duration T of 
the dataset, the mode amplitude A and the background, 
power density B (including stellar and white noise) at the 
frequency of the mode. In the case of both Sun-like and 
red giant oscillations, mode lifetimes are of the order of a 
few days (Fletcher et al. 2006), much shorter than a typical 
CoRoT run (21 or 150 days). We therefore use the formula 
given by Libbrecht (1992) for the case where r <C T to 
estimate the frequency uncertainty o v as: 



oi = /(/?) 



r 



(2) 



where /3 = B/A is the inverse signal-to-noise ratio, and 



/(/3) = (1 + /?) 1/2 [(i + ^Va + zJi/a 



As is common practice in asteroseismology (see e.g. 
Kallinger et al. 2009 and references therein), we model B 
as a sum of broken power-laws of the form 



+ (27rbiv) c * 



(3) 



where each term in the sum corresponds to the stellar vari- 
ability associated with a particular type of surface struc- 
ture (typically active regions and granulation, with in some 
cases a super-granulation component), aj is the amplitude, 
bi the characteristic time-scale, and Ci the slope of term i, 
and d represents the white noise. 

As a proxy for the background signal produced by a 
typical star in which one might search for Sun-like oscilla- 
tions, we fit a 3-component background (following Eq.[3]) to 
total solar irradiance data obtained by the VIRGO/PM06 
radiometer on the SoHO satellite during a 5-month period 
close to the last solar activity cycle maximum (for details 
see Aigrain et al. 2004). We then added white noise at the 
level expected for a bright (V = 6) star in the CoRoT 
asteroseismology (AS) channel (all white noise values from 



Auvergne et al. 2009). The resulting background power den- 
sity spectrum is shown as the thick solid line in the top 
panel of Fig. [1] over the frequency range where oscillations 
are expected. We also constructed a less favourable case by 
scaling up the solar background by a factor of 10 (this is 
not unreasonable, as the Sun is generally thought to be a 
relatively quiet star) and adding white noise corresponding 
to a V = 8 star observed by CoRoT (thick dashed line in 
Fig. [1]) . Finally we also constructed a red-giant case by vi- 
sually estimating the power density spectrum parameters 
for such an object from Fig. 1 of Kallinger et al. (2009) and 
adding white noise corresponding to a V = 12 star observed 
by CoRoT (thick dash-dot line in Fig. [TJ). The adopted 
background power density parameters are summarised in 
Table [TJ 

We then computed frequency uncertainties as a function 
of mode frequency, using Eq. (|2|), for mode amplitudes of 
A = lOppm (sun-like oscillation in AS channel) or A = 
1000 ppm (red giant oscillations in PF channel). The results 
are shown in the bottom panel of Fig. [TJ In both cases, we 
assumed a mode lifetime of T — 5 days and a run duration 
of T = 150 days (CoRoT long run). 



150 d run, mode lifetime = 5 d 



Giant, PF channel, V=12 



10 x active Sun, AS channel, V=8 _ 
active Sun, AS channel, V=6 



mode amplitude 10 ppm 
mode amplitude 1000 ppn 



100 1000 
Frequency (^Hz) 



Fig. 1. Top panel: Typical background power spectrum 
density B{v) (see text for details) for a bright Sun- like star 
in the CoRoT AS channel (solid curve), a fainter and more 
active star in the same channel (dashed curve), and a red 
giant at the bright end of the PF channel (dash-dot curve). 
The thinner horizontal lines show the white noise contribu- 
tion in each case. Bottom panel: resulting frequency uncer- 
tainties for modes with a lifetime of 5 days observed for 150 
days, with amplitudes 10 ppm (AS channel) and 1000 ppm 
(PF channel). Also shown are the minimum uncertainty 
arising from the finite mode lifetime (horizontal solid grey 
line) and the minimum uncertainty arising from the finite 
run duration for an infinite mode lifetime (Nyquist sam- 
pling, horizontal dashed grey line). 
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3. Analysis of theoretical oscillation spectra 

We first generate a coarse grid of equilibrium stellar mod- 
els with various values of mass (1.35 < -M- < 1.45, 4¥- = 

0.01), luminosity (0.55 < log(j^) < 0.65, S(log(^)) = 
0.01) and effective temperature {6300K <T eff < 6500K, 
5T e ff — 25K). A refined grid is then generated around 
the stellar parameters which predict close pulsation spec- 
tra, according to the results obtained from the coarse grid. 
On this refined grid (S(log(^-)) = 0.005 , 5T eff = 10K, 

but constant mass because of time computation), we also 
vary the mixing-length parameter (1.6 < oimlt < 1-8, 
8 oimlt = 0.05), the metallicity and helium content Ya e 
(0.26 < Y He < 0.28, SY He = 0.005) of the equilibrium 
model. For this preliminary study, the metallicity Z is var- 
ied between Z = 4.10 -3 and Z = 0.02. In a forthcoming 
study, we plan to include a finer grid for the stellar mass 
and models with over-solar metallicities. 

We use two methods for the comparison between pul- 
sation spectra computed from two different sets of stellar 
parameters: 

— the maximum difference of the frequencies dv of the 
spectra for given values of the radial order n and of the 
degree I of the mode, 

— the maximum difference of the large separation Av = 
v n ,i — v n -i.i or small separation 5vq2 = ^«,o — ^n-1,2, 
(Kjeldsen et al. 2008). 

A comparison of these two methods is given below. 

3.1. Effect of the mixing length parameter olmlt 
and of the metallicity Z 

We examine the effect of the mixing length parameter 
olmlt on the spectra in the following way : we fix M, L, 
T e ff, Yff e and Z and compute, for a given mode (n, I), the 
difference dv between the frequency for a value of cvmlt 
and the frequency for a reference value of oimlt (in this 
case olmlt, re f = 1-7). The results are shown in Figure [2] 
We performed a similar exercise to estimate the effect of 
the metallicity Z. The results are displayed in Fig. [3] using 
a reference value for the metallicity Z re f=0.01. Our prelim- 
inary tests show that a variation of oimlt or of Z, within 
reasonable values, can yield a variation of the frequency of 
a given mode of several fiHz. These effects need to be in- 
vestigated in a more systematic study, which we plan for a 
forthcoming work. 

3.2. Comparison between abolute values of frequencies 

For a given relative mass difference 5M/M, the minimum 
value Sv m in of 5v is determined in the coarse grid. The 
value of 5v m in depends on the number of considered modes. 
We vary this number from 4 (one mode for each value of 
< I < 3) to 60 (15 modes for each value of < I < 3). 
The variation of Sv m i n with 8M/M, for different numbers of 
considered modes, is shown in figure[4] The same procedure 
is applied for a given luminosity difference SL/ L or a given 
radius difference SR/R. The results are shown in figure O 
and figure [6] respectively. 

As expected, the larger the mass/luminosity /radius dif- 
ference, the larger the spectra difference. Similarly, the 
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Nb of node 

Fig. 2. Frequency difference dv (/iHz) as a function of 
the radial order n, for different values of oimlt compared 
to the reference value oimlt — 1.7. The solid red/dashed 
blue/dotted green/dash-dotted black curves correspond to 
degree I of the mode equal to 0/1/2/3 respectively. The hor- 
izontal line is a guide line which corresponds to a frequency 
difference dv=0. 



15 r 




Nb of node 

Fig. 3. Frequency difference dv (//Hz) as a function of the 
radial order n, for different values of metallicity Z com- 
pared to the reference value Z = 0.01. The solid red/dashed 
blue/dotted green/dash-dotted black curves correspond to 
degree I of the mode equal to 0/1/2/3 respectively. The hor- 
izontal line is a guide line which corresponds to a frequency 
difference dv=Q. 
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4^ 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 
5m/m 

Fig. 4. Spectra difference Sv (//Hz) as a function of the rel- 
ative mass difference SM/M for the coarse grid of stellar 
models (see £)3.2[) . The different curves correspond to dif- 
ferent numbers of considered modes, as indicated near the 
curves. 



larger the number of considered modes, the larger the spec- 
tra difference. These trends are the same on the refined grid, 
but the values for Sv m i n are much smaller. If we focus on 
the mass difference, for instance for 40 modes, and SM/M 
— 0.06, SiSmin decreases from 3.4 /tHz on the coarse grid to 
0.6 /iHz on the fine grid, as illustrated in Fig. [7l The effect 
is similar for the luminosity difference (for 40 modes, and 
SL/L — 0.06, Sumin decreases from 2.4 /iHz on the coarse 
grid to 0.5 /iHz on the fine grid) or the radius difference 
(for 40 modes, and SR/R = 0.03, 8v m in decreases from 6.4 
/iHz on the coarse grid to 1.0 /(Hz on the fine grid). Since 
we cannot freely refine the grid because of computing time, 
the values obtained for 8v m i n on the refined grid may be 
overestimated : we cannot exclude that a finer grid gives 
smaller values of Sv m i n . 

3.3. Large and small separation 

Many asteroseismological studies do not directly compare 
the pulsation frequencies of two stars but rather compare 
the large and small separations. For this second method, 
we define the difference 5 s i s between two spectra as the 
maximum value of the large separation or of the small sep- 
aration (the highest value is retained). The results for the 
coarse grid, for the mass difference 8M/M with 60 con- 
sidered modes, are shown in Fig. [5] S s i s is about a factor 
of two lower than Sv. For the luminosity difference, S s i s is 
also about a factor of two smaller than Sv. For the radius 
difference, S s i s is about a factor of three lower than Sv. 

As in paragraph 13. 2[ the same behaviour is found when 
using the fine grid for a given SM/M, SL/L or SR/R. The 
results corresponding to a given SM/M are shown in Fig. 

m 
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o — 1 1 1 ' 1 — 1 1 1 1 1 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
5L/L 

Fig. 5. Spectra difference Sv (/iHz) as a function of the 
relative luminosity difference SL/L for the coarse grid of 
stellar models (see Ej3.2[) . The different curves correspond 
to different numbers of considered modes, as indicated near 
the curves. 

We also examined a third method, where the difference 
< 8 s is > between two spectra is defined as the difference 
between the mean large separation or the mean small sep- 
aration (the highest value is retained) . We did not find this 
method relevant in the present context (< 5 s i s > is found 
to be negligible for small values of SM/M, SL/L or SR/R 
in the grids studied). 

4. Discussion and conclusion 

The results obtained in <J3] enable us to derive the accuracy 
of the stellar parameters mass M, luminosity L and radius 
R determined from asteroseismological data. The accuracy 
for one parameter is obtained letting all the other parame- 
ters remain free. According to £12.31 we adopt the values of 
0.2 /iHz and 0.3 /iHz respectively for the precision of fre- 
quency determination based on the CoRoT performances. 
The smaller value is for a red giant in the PF channel and 
the larger value corresponds to the accuracy obtained for a 
Sun-like star in the AS channel (see Fig. [TJ). 

The results are given in Tables [HO and Q] for an accuracy 
of the frequency determination of 0.2 /iHz, and in tables [5] 
[6] and [7] for an accuracy of the frequency determination of 
0.3 /tHz. The most stringent constraint is provided by the 
method using Sv. For Sun-like stars as observed in the AS 
channel, our analysis suggests for instance an uncertainty 
for the mass as large as 5% if only 20 modes are considered 
(for the closest models, effective temperatures are compati- 
ble within 60K, and radii and luminosities differ by ~ 2%). 
With the method based on S s i s , an uncertainty of 5% for the 
mass remains for up to 50 considered modes (for the closest 
models, effective temperatures are compatible within 60K, 
radii differ by 3% and luminosities are equal). 
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0.005 



0.01 



0.015 

5R/R 



0.02 
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0.03 




30 40 
Nb of mode 



Fig. 6. Spectra difference 8v (//Hz) as a function of the rel- 
ative radius difference SR/R for the coarse grid of stellar 
models (see £)3.2[) . The different curves correspond to dif- 
ferent numbers of considered modes, as indicated near the 
curves. 

Table 2. Accuracy SM/M on mass determination, for a fre- 
quency accuracy of 0.2 /iHz, obtained from different meth- 
ods of evaluating spectra difference and for different values 
of considered modes. 



Nf of modes 


20 


40 


60 


SM /M with Sv 


0.02 


0.01 


<0.01 


SM/M with 5 sis 


> 0.07 


0.02 


0.01 



Table 3. Accuracy 8L/L on luminosity determination, for 
a frequency accuracy of 0.2 /iHz, obtained from different 
methods of evaluating spectra difference and for different 
values of considered modes. 



Nf of modes 


20 


40 


60 


SLj L with Sv 


0.06 


0.04 


0.02 


SL/L with 8 s is 


0.1 


0.05 


0.05 



Table 4. Accuracy 8R/R on radius determination, for a fre- 
quency accuracy of 0.2 /iHz, obtained from different meth- 
ods of evaluating spectra difference and for different values 
of considered modes. 



Nf of modes 


20 


40 


60 


8R/R with Sv 


<0.01 


<0.01 


<0.01 


5R/R with 8 s is 


0.02 


0.01 


<0.01 



Our work is a first step toward a more systematic study, 
as it only takes into account variations in the basic input 
physics of stellar models such as the mixing length param- 
eter, the helium abundance and to a certain extent the 
metallicity. A more systematic analysis should also explore 
the effects of different convection models, of core overshoot- 



Fig. 7. Spectra difference 8v (see §3.2j) as a function of 
the number of modes for the fine grid of stellar models. 
The different curves correspond to different values of the 
relative mass difference SM/M. 



Table 5. 

/iHz 



Same as Table [5] for a frequency accuracy of 0.3 



Nf of modes 


20 


40 


60 


SM/M with Sv 


0.05 


0.02 


0.02 


SM /M with 8 sis 


> 0.07 


0.06 


0.04 



Table 6. Same as Table [3] for a frequency accuracy of 0.3 
/tHz 



Nf of modes 


20 


40 


60 


SL/L with Sv 


0.09 


0.05 


0.03 


SL/L with 8 sis 


>0.1 


>0.1 


>0.1 



Table 7. Same as Table [4] for a frequency accuracy of 0.3 
/tHz 



Nf of modes 


20 


40 


60 


8R/R with Sv 


0.01 


<0.01 


<0.01 


5R/R with Ssis 


0.025 


0.015 


0.01 



ing and of rotation. These uncertainties in current stellar 
evolution theory will certainly increase the degeneracy of 
the pulsation spectra (i.e the limit on the mass, luminos- 
ity or radius mentioned above may be significantly larger). 
Moreover, we have not explored the effect of uncertainties in 
the stability analysis calculation, such as the sensitivity to 
boundary conditions. The large and small separations are 
expected to be less sensitive than absolute frequencies to 
boundary conditions used in the computation. Obviously, 
much additional works need to be done. Our work, how- 
ever, provides an idea of the level of accuracy required by 
asteroseismological studies in order to improve exoplanet- 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 

5m /m 

Fig. 8. Comparison between two methods measuring spec- 
tra differences for the coarse grid of models and with 60 
considered modes. Red/solid line : frequency difference 
( §3.2j) . blue/dashed line : large or small separation differ- 
ence ( %T3|) . 



0.8 




0.1 - 

Q l 1 1 1 1 1 

10 20 30 40 50 60 

Number of modes 

Fig. 9. Comparison between two methods measuring spec- 
tra differences for the fine grid of models and for SM/M 
= 0.06. Red/solid line : frequency difference ( fl3.2[) . 
blue/dashed line : large or small separation difference 

(833). 



host star parameters, crucial information for future projects 
such as PLATO. 
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Case 


V 


A 


d 


ai 


bi 


Cl 


a-2 


b 2 


C2 


a.3 


b 3 


C3 


active Sun 


6 


10 


0.65 


2.0 ■ 10 4 


2.0 • 10 b 


4.35 


5.51 


8205 


3.5 


0.19 


341 


1.4 


10 x active Sun 


8 


10 


1.95 


2.0 ■ 10 5 


2.0 • 10 6 


4.35 


55.1 


8205 


3.5 


1.9 


341 


1.4 


red giant 


12 


1000 


6.8 


3.5 ■ 10 5 


1.2 • 10 5 


4.0 


55.1 


6367 


3 









Table 1. Parameters of the background models used. The photon noise level d is given in ppm/V^Hz, while the 
background amplitudes ai are given in ppm 2 /^Hz, and the time-scales bi in 10 6 s. Three components were used for the 
Sun-like cases but 2 only for the red giant case. 



